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We examine the stop production scenario for the anomalous excess of high Q 
events observed at HERA, by taking into account the constraints coming from 
the leptoquark search at Tevatron and the atomic parity violation experiments. 
This scenario is shown to survive persistently even under the severe constraints 
from the other experiments. In the analysis, the branching ratios of the 

decay modes tj — > ed (J = 1, 2) are found to be useful parameters, representing 
a number of unknown parameters in the SUSY models. By expressing the stop 
contribution to the cross section of deep inelastic ep scattering a(ep — > eX) 
in terms of B^ and B^ 2 , the allowed region is successfully identified in the 
parameter space of B^ and B^ 2 . 
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An anomalous excess of high Q 2 events in the deep inelastic scattering e + p — > e + X at 
HERA has been reported by HI and ZEUS collaborations^, |2|, [3]. The excess is observed 
for high momentum transfer of positron, Q 2 > 15000GeV 2 , and for relatively high x region, 
x > 0.25, where x is the Bjorken scaling variable. The measured cross section for the 
excess events is 0.71pb against the standard model(SM) expectation of 0.49pb. In addition, 
the excess seems to be broadly distributed around 200 ~ 250GeV for the invariant mass 
distribution of the positron-jet system. In the following, let M be the invariant mass of 
positron-jet system represented by M = ^/xs, where the center of mass energy of the e + p 
system = 300GeV at HERA. 

In order to explain the excess of high Q 2 events, several scenarios have been proposed j|, 
HI H> & 0- Among them the most plausible scenarios are the squark production ones in 
which the excess and the M distribution of these high Q 2 events are caused by the s-channel 
production of the squark following the positron-quark collision. Although the high Q 2 events 
can be enhanced by a s-channel resonance, the shape of the M distribution by the s-channel 
production of a single squark is too sharp to explain the data. Therefore, it is difficult to 
explain the broadness of the M distribution by a single squark production.^ Fortunately, in 
the squark production scenarios, two almost degenerate peaks can be arisen around M = 
200 ~ 250GeV, and therefore the moderate M distribution may be explained due to the poor 
statistics of the current data at HERA experiments^]. Furthermore, although the squark 
production scenarios are possible to explain the features of the excess of high Q 2 events at 
HERA, we must take into account the constraints coming from both high and low energy 
experiments. 

In this article, we consider the squark production scenarios in the supersymmetric(SUSY) 
model with i?-parity breaking^) interaction and examine whether the squark production 
scenario satisfies the constraints from the current experiments. We show that the squark 
production scenario is consistent with the constraints from the leptoquark search at Tevatron 
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1 The discussion on the single squark production scenario can be translated to that for the scalar lepto- 
quark production scenario with suitable substitution of coupling and mass. 
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and the measurements of atomic parity violation. Also it is pointed out that the useful pa- 
rameter is the branching ratio Bq of the decay mode q — > eg, when we discuss the consistency 
of the squark production scenarios with the leptoquark search at Tevatron. The constraint 
on the plane of the squark mixing parameter 6g and the Iji coupling A' is also discussed. 

The existence of the lepton-squark-quark(7-g-g) couplings is essential to explain the excess 
and the M distribution of high Q 2 events. The l-q-q couplings are induced by Iji interaction 
that results from the trilinear ]fl superpotential.f] The general form of the $ superpotential 
is given by 



Wft = + K 1k UQ 3 bt + \" jk u?Dw 



where L and Q are the lepton and quark superfields in SU (2) doublet representation, whereas 

E c , U c and D c are superfields for lepton, up- and down-type quarks in singlet representation, 
respectively; and i, j and k denote the generation indices. The Iji interaction Lagrangian 
relevant to the high Q 2 excess events is derived from Wft with i — k — 1 as follows: 

C = X XjX {u jL dP L e - d R e c P L Uj) + h.c, (2) 

where Pl is the left handed chiral projection operator. In general, the squarksf^, qn) are 
represented as the mixed states between the mass eigenstates (q±, q 2 ) of squark with a mixing 
parameter 8g as follows: 

q L = qi cos 9$ - q 2 sin 0%, qR = <?i sin 6% + q 2 cos Oq. (3) 

As mentioned above, two almost degenerate peaks around 200 ~ 250GeV are required to 
explain the broadness of M distribution. The squark mixing is necessary to produce two 
peaks in the M distribution by ]/L interaction. 

Hereafter, as a representative of the squark production scenarios, we consider the stop 
production scenario, in which case j = 3, and q should be identified as the stop(t) f\. From 
the combined data of both HI and ZEUS, the excess events with Q 2 > 15000GeV cluster 
around M ~ 200GeV and M ~ 230GeV.Q The location of the peaks in the M distribution 
should be identified as the mass of the two stops. Hereafter we take = 200GeV and 
m~ t2 = 230GeV. 

There exist many constraints on the absolute value of ]ft coupling constants or on the 



products of themfTT|, |T2]j . Especially, the strict constraints come from the proton lifetime p"3|| 



neutrinoless 2(3 decay P^j , neutrino mass bound(m- )\T§\, K -K°(D° -D and B°-B°) mixing 



and both rare and forbidden decays of mesons Pq| , fi-e conversion ||17|| and atomic parity 
violation(APV)lT§. 

The non-zero value of A' 131 is necessary to explain the excess events at HERA by the stop 
production scenario. Hereafter, we suppose only A' 131 to be non-zero and other Iji coupling 
constants are to be zero. As a result we can relax all the above constraints except the one 
from APV. 



2 On the other hand, leptoquark is a hypothetical particle that couples to both lepton and quark at the 
same vertex. 

3 Scalar charm quark(c) production scenario can be considered when we take j = 2 and q = c[^|, ^, 

The stop production scenario can be paraphrased into the scharm production scenario by suitable substitution 
of mixing angle, coupling and mass. Scalar up quark production scenario is, however, impossible to explain 
the event excess at HERA. It is because the coupling constant of X' nl is already restricted within a very 
small value by neutrinoless double beta decayQ for the scalar up quark mass around 200GeV and gluino 
mass smaller than a few TeV, and as a result its contribution to the excess of high Q 2 events at HERA is 
very small. 

4 Since the statistics are not high enough to use HI and ZEUS data separately, we use the combined data 
of both HI and ZEUS. One may doubt whether the events clustered around ~ 230GeV can be identified as 
the peak, though it may be possible that the events clustered around ~ 200GeV to be identified as the peak. 
In this article, we will consider the case that the two peaks exist. 
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The most precise measurements of the weak charge Q e ^ p in the APV experiments have 
been obtained from 133 Cs. The weak charge is measured to be Q e ^ = — 72.11±0.27±0.89[|T9"|. 
The SM prediction including the radiative correction [pOfl is given as Qy/ 1 = —73.11 ± 0.05 
leading to AQ W = Q e ^ p - = 1.00 ± 0.930. The bound on AQ W with 95% C.L. is 
obtained as follows: 



- 0.82 < AQ W < 2.8. 
The contribution of the stop to the weak charge JT^] reads 



AQ 
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where N = 78, Z = 55 for 133 Cs. In figure [l], the upper bound on A' 131 is shown as a function 
of cos% We note from figure [l] that A' 131 < 0.07 is allowed in the whole region of &i when 
m h = 200GeV and m h = 230GeV.Q 

The formula for the cross section a{ep — > eX) with double stop production can be derived 
from the straightforward extension of the formula with single stop production in the^£-SUSY 
model[10]. Besides the SM parameters, a(ep — > eX) depends on six variables 



<j{ep — > eX) 



v , m t 2 , Oh A 'l31 > T h > T h ) > 



(6) 



where are the total decay widths of tj. When it is kinematically possible that the tj 
decay modes include SUSY particles, ry may depend on a number of SUSY parameters.^ 

Once we treat the branching ratios as independent parameters, however, the contribution to 
a(ep — > eX) from the parameters of SUSY model, except for those of the stop sector, can 
be embedded into the two finite parameters, < Br(tj — > ed) < 1. 
The total decay widths can be expressed as follows: 



F(tj — > ed) 
Br{tj — > ed) 



(J = l,2), 



(7) 



where T(tj — > ed) are the partial decay widths of the stop to the ed mode, and Br(tj — > ed) 
are the corresponding branching ratios. The formulae of T(tj — > ed) are 



r(t! -> ed) 



I Ksi 1 2 cos 2 flf 
16tt 
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Hereafter we will use shortened notation Br 
rewritten as a function of the six parameters, 

o{ep — * eX) 



r(t 2 -> erf) 



|A' 131 | 2 sin 2 % 



rn 
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erf). Then o{ep 



<y K , m t - 2 , Ai 31 , B h , s t - 2 ; 



eX) can be 



(9) 



This choice of parameters is convenient when we discuss the feature of event excess at HERA 
almost independent of the detailed structure of SUSY model. 

The branching ratios are restricted directly by the first generation scalar leptoquark(LQi) 
search at Tevatron experiments. The upper bound on the a(pp — > LQ\LQ 1 X)Bj j Q i is 

obtained by the Tevatron exp eriment s ||2 2[| , where a(pp — > LQiLQiX) is the total cross 
section of the LQi pair production and 1?2qi is the branching ratio of the LQi — > eq decay 

5 The constraints on Ifl interaction may be relaxed by the introduction of other new physics sources, for 
example, the contact interactions . 

6 For example, when the decay mode t — * x°Q is kinematically allowed, the total decay width depends 
on gaugino mass parameters (M^ (i = 1,2)), Higgsino mass parameter(/i) and the ratio of the vacuum 
expectation value of two Higgs bosons(tan/3 = = ^), in addition to m t ~, 0x and A' 131 . 
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mode. The upper bound can be interpreted as that of a(pp —> ttX)B 2 , because one expects 
that the signal selection cuts for LQi can be adopted for the stop with^? interaction ([|). The 

theoretical prediction of a(pp — > ttX) has been estimated including SUSY QCD effects at the 
next-to-leading order Consequently, we can obtain the upper bounds on the branching 
ratio Bf. The upper bound are listed in Table [I]. 

Figure ^| shows the contour plots of a(e + p — ► e+X) in the B$ versus B^ 2 plane. Although 
the branching ratios are constrained as listed in Table |1] by the first generation leptoquark 
search at Tevatron, there remains a region in the Bi x -B~ t2 plane where the excess of high Q 2 
events at HERA can be explained by the stop production scenario. We also find out from 
figure [| that the contours of the total cross section a(e + p — *> e+X) for the same value of 9i 
seem to be parallel. The reason why the contours are parallel in figure ^| can be explained 
as follows. 

The approximate formula of da(e + p — > e + X)/dxdQ 2 is obtained from its exact formula 
by using the narrow width approximation; 

da(e+p -> e+X) da(e + p -> e+X) 



dxdQ 2 



dxdQ 2 
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where the first term on the right hand side, denoted by da ^ x 9 d 

t-channel exchange of 7, Z and tj (j = 1,2), whereas the second and third terms come from 
the s-channel exchange of t\ and t2, respectively. The parton distribution function for the 
down quark inside a proton p5| is denoted as d(x, Q 2 )- Since the amplitude via the t-channel 



exchange of the stop, as compared with that by the SM process for m^ > 200GeV and 

1I0D is approximately 
formula for the total 



A' 131 < 0.1, is negligibly small, the first term on the right hand side in 
equal to the contribution from the SM. Therefore, the approximate 
cross section a(e + p 
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e r X) can be written as follows: 
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where /(m? / 's,Q 2 mir ) are defined by 



I{m%/s, Q 2 min ) = f dx I™ dQ 2 d(x, Q 2 ) 5(sx - m? ) (j = 1, 2) 



Xmin = Qminl s i an d Qmin * s t ne selection cut for the high Q 2 events. Equation(|TT|) shows 
that cr(e+j9 — > e+X) depends on B^. (j = 1, 2) linearly. This explains the feature of figure [2] 
that the contours of the total cross section are almost parallel in the Bi -Bi plane. 
After integrating over x, we have 



7(m?/s, Q 



2 ■ ; 

mm ) 
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min 



dQ 2 d{x, Q 2 
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= 1,2). (13) 



(if m? < QLn) 
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Several explicit values of I(m~ /s,Q^ n ) are given in Table |[ From eq.flTTD and Table 0, we 
can easily calculate a(e + p — > e + X) for arbitrary values of A' 131 , % 1?^ and I?j 2 by substituting 
the two values for j — 1, 2 of I(m 2 Js, Q^m) with the same Q^m- I* 1 figure |3|, we show the 

contour plots of a(e + p —>■ e + X) derived by using the approximate formula(dashed lines) and 
the exact formula(solid lines). 

Next, we discuss whether the constraint on the 9i and A' 131 can be obtained from the data 
of HI and ZEUS for excess events. Combining data of HI and ZEUS for Q 2 > 15000GeV 2 , 
we can find 11 events in the first bin 187.5GeV< M < 212.5GeV and 5 events in the 
second bin 222.5GeV< M < 247.5GeV with luminosity of L = 57.2 ph' 1 ^ against the SM 
expectation of 5.9 events in the first 187.5GeV< M < 212.5GeV and 2.1 events in the second 
222.5GeV< M < 247.5GeV. We now require the constraint on the number of events in each 
bin at la level as follows: 

7.7 < iV| 187 .5 G eV<M<212.5GeV < 14.3, 2.8 < N\ 222.5GeV<M<247.5GeV < 7.2. (14) 

Figures|](a) and |](b) show the contour plots of the expected number of events for several sets 
of Bi x and B^ 2 , along with the upper bound on A' 131 obtained by APV measurements.^ In 
figures^j(a) and f|(b), the constraint on the cos^-A' 131 plane becomes more strict for smaller 
values of B^ . It is because the total decay width of t becomes larger for smaller values of B^ , 

and as a result the total cross section of a(e + p — > e + X) decreases with decreasing Bi for a 
fixed value of A' 131 . Therefore, to explain the excess of events, we must take the larger values 
of A' 131 corresponding to the smaller values of B% (j = 1,2). The value of A' 131 , however, 

cannot be taken too large for the small B^, since A' 131 larger than 0.07-0.08 is excluded by 
APV. 

To conclude our discussion, we show that the stop production scenario for the excess 
of high Q 2 events at HERA is consistent with the current experiments, even though the 
leptoquark search at Tevatron and APV measurements strictly constrain the stop production 
scenario. 

Currently, the statistics at HERA experiments are not enough to constrain the cos 6^-A' 131 
parameter space strictly. In the future high luminosity run at HERA, we expect that the 
number of events in each bin will be restricted within a narrower region than the present 



one in eq. (|14D , so that a more strict constraint on the cos#£-A' 131 parameter space will be 
obtained in the near future. HERA upgrade is planned with L w 7.4 x 10~ 5 pb~ 1 s~ 1 starting 



in the year 2000 [26] Furthermore, if a large amount of luminosity is accumulated at the 



future Tevatron experiment such as TEV33(L = 210.62pb 1 / week)|27] by the end of year 
2006, then we may expect that the branching ratios B^. (j = 1, 2) will be restricted within 

the small values, namely B^ < 0.2 and B^ 2 < 0.33. Therefore the stop production scenario 
will be strictly constrained by TEV33. Nevertheless, the figures E] and |] show that the stop 
production scenario can be compatible with the constraints from both the leptoquark search 
at TEV33 and the measurements of APV, even if the strict constraint on the B^ (j = 1, 2) 
is obtained at TEV33. 

The authors thanks M. Ahmady for reading the manuscript and helpful comments. One 
of the authors (J. K.) thanks T. Kon and F. Shibata for helpful discussion and comments. 
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Figure captions 



Figure [T] The upper bound on the value of A' 131 from the measurements of the atomic 
parity violation(APV) is shown as a function of cos% We take (m^, m^ 2 ) = (200GeV, 
230GeV). 

Figure [2| Contour plots of the cross section of a(e + p — ► e + X) for Q 2 > 15000GeV 2 are 
shown in the B^ versus B^ 2 plane; the dashed line is for cos^ = 0.3, the solid line 
for cos#£ = 0.6 and the dotted line for cos^- = 0.9. We take A' 131 = 0.07 and y/s = 
300GeV. Each value written beside the curve is corresponds to the measured value, 
o(e + p -> e+X) = 0.71±g;i|pb, for Q 2 > 15000GeV 2 . 

Figure [3| Contour plots of the cross section of a(e + p — > e + X) for Q 2 > 15000GeV 2 and 
cos^j- = 0.6 are shown in the B% versus Bi 2 plane; the dashed and solid lines are 
corresponds to the results obtained by the approximate formula and the exact formula, 
respectively. Other input parameters are the same as in figure 0. 

Figure H(a) (b) Contour plots of the expected number of events in each bin; 187.5GeV< 
iW<212.5GeV and 222.5GeV< M <247.5GeV; are shown in the cos0 t -A' 131 plane, 
where \fs = 300GeV. As for the branching ratio of tj (j = 1,2), we take (a):(Sf , 
B^)=(0.65, 1.0) and {b):(B h , B h ) = (0.2, 0.33). The integrated luminosity is taken to 
be 57.2pb _1 . The curve of the upper limit on A' 131 by APV are also superposed. 
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Tables 





experimental 


experimental 


theoretical 


mass 


upper bound on 


upper bound on 


prediction |23 


m t 


Br(t -> ed) 


a(pp -> ttX) {Br (t -> erf)} 2 


-> iiX) theor , 


200 GeV 


< 0.65 


< 0.076 pb 


0.18 pb 


230 GeV 


< 1.0(no bound) 


< 0.067 pb 


0.065 pb 



Table 1: The upper bound on Br(i — > ed) and cr(pp — > iiX){Sr(t — > ed)} 2 
from the first generation leptoquark search at Tevatronp2[. 



mi 


QLn = lOOOOGeV 2 


Q 2 mtn = 15000GeV 2 


180 GeV 


7.061 10-2 


5.438 10-2 


190 GeV 


5.561 10-2 


4.456 10-2 


200 GeV 


4.160 10-2 


3.434 10-2 


210 GeV 


2.937 10-2 


2.481 10-2 


220 GeV 


1.937 10-2 


1.668 10-2 


230 GeV 


1.177 10~2 


1.028 10~2 


240 GeV 


6.437 10-3 


5.694 10-3 



Table 2: Values of I(m,f/s, Qmin) f° r several values of the stop mass(m 4 ~) and 
Qmin- Definition of I(m-J s,Q 2 min ) is given by eqs.(H) and @). 



9 



^0.7 



009 z_ excluded by APV (95% CX.) 

0.08 ^ 

0.07 ^ 

0.06 ^ 

0.05 ^ 

0.04 ^ 

0.03 ^ 

0.02 ^ 

0.07 - 11^2 = 230 GeV 



m- tl = 200 GeV 



0.7 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cos§~ 



Figure 1: 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Br(tj-> ed) 



Figure 2: 
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Br(tj-> ed) 



Figure 3: 
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Figure 4: (a) 
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Figure 4: (b) 
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